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Abstract 


A  variational  approach  is  applied  to  collisions  of  a  particle 
having  a  mass  tn  with  a  harmonic  oscillator  having  a  mass  M  in  one 
dimension.  We  have  used  trial  functions  which  become  exact  in 
the  limit  of  either  m/M>>l  or  m/M«l  .  We  have  derived  exact 
expressions  in  the  limit  m/M  —  0  for  inelastic  as  well  as  elastic 
collisions  to  leading  order  in  the  parameter  m/M  for  any 
impenetrable  interaction  potential. 
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I .  Introduction 


Scattering  of  atoms  from  surfaces  has  received  considerable 
attention  in  recent  years.1  Several  theoretical  methods  have  been 
developed  to  deal  with  this  problem:  The  stochastic  classical 
trajectory  method,2  the  eikonal  approximation,3  the  trajectory 
approximation,4  and  the  semiclass ica 1  gaussian  basis  set  method.5 
So  far,  however,  fully  quantum  mechanical  treatments  are  limited 
to  theories"  based  on  (first-order)  distorted  wave  Born 
approximation  (DWBA) .  In  this  paper,  we  propose  another  quantum 
mechanical  treatment  by  employing  a  time-independent  variational 
approach . 9 

Variational  methods  in  scattering  problems  can  be  classified 
into  three  groups:  The  first  one  is  based  on  the  standard 
variational  principles  of  Hulthen10  and  Kohn.11  The  second  one  is 
based  on  the  Schwinger  variational  principle.12  The  last  one  is 
the  combination  of  these  two.  first  formulated  by  Takatsuka  and 
Mckoy13  for  the  phase  shift  and  then  by  Gross  and  Runge14  for  the 
T-matrix.  The  advantage  of  this  approach  is  that  we  can  evaluate 
the  T-matrix  -as  easy  as  we  can  in  the  first  approach  without 
imposing  the  standard  scattering  boundary  conditions  on  trial 
functions.  For  this  reason,  we  will  take  the  third  approach. 

In  the  present  work,  we  consider  a  quantum-mechanical  system, 
in  one  dimension,  of  a  particle  interacting  with  a  fixed  harmonic 
oscillator.  This  is  a  paradigmatic  model  and  provides  a  ground 
for  testing  any  approximate  method  to  treat  the  problem  of 
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atom-surface  scatterings.  The  Hamiltonian  of  the  system  is 


d ^  h~  d~  ,  M  2  2.1// 

2m  ax2  2>l  3y-  2  “  y  V(X  y)  ' 


where  m  and  x  are.  respect i vely ,  the  mass  and  coordinate  of  the 
incident  particle  "incidon'),  M,  a,  and  y  are  the  mass, 
frequency,  and  coordinate  of  the  oscillator  r espect i vel y ,  and 
V ! x-y)  is  the  mutual  interaction  which  we  assume  impenetrable. 
When  we  measure  energies  and  lengths  in  units  of  hcj  and  ( h/A/co ) ,/2 
respect i vely ,  the  Hamiltonian  (1)  can  be  rewritten  as 


u  _ _ 1 

2pax2 


i  +^y9  +  V'(x-y) 
^  dy~  ^ 


where  p  is  the  mass  ratio  m/M. 


with  the  numerical  results  was  obtained  by  a  simple  semiclassical 


S-matrix  treatment.  The  method  of  Heller5  could  also  give  as 
accurate  results  as  Miller's  method  could. 

The  criterion  for  the  validity  of  semiclassical  calculation 
is  that,  for  given  initial  oscillator  and  incident  energies,  the 
possible  classical  energy  changes,  AEct  (which  depend  on  the  angle 
parameter  of  the  oscillator),  span  several  oscillator  excitation 
quanta : 


An  -  AEd/Ku  »  1  .  (2) 

This  is  an  expression  of  the  correspondence  principle.  As  is  well 
known,  suitable  semiclassical  calculations  may  give  rather 
accurate  results  even  when  the  condition  (2)  is  not  satisfied20. 

For  given  initial  energies  of  the  oscillator  and  incidon. 
the  semiclassical  condition  (2)  is  violated  in  two  regimes  when. 


as  is  well 

known , 

classical  energy 

transf  er 

is  smal 1 :  v  I ) 

m»M 

and 

(II) 

m«M . 

The  first  regime 

is 

similar  to 

the 

Born 

-Oppenheimer 

adiabatic  regime 

for 

molecules ,  with 

t  he 

oscillator  and  incidon  corresponding  respectively  to  the  molecular 
electrons  and  nuclei.  The  second  regime  may  be  called 
anti-adiabatic.  We  will  treat  primarily  these  two  extreme 
situations  in  this  paper. 

In  both  regimes  the  collision  is  predominantly  elastic.  The 
inelastic  collisions  depend  on  the  details  of  the  potential 
V  ( x-y)  i-i  the  first  regime,  but  they  do  not  in  the  second  one 
Thus  we  can  derive  the  exact  expressions  for  inelastic  collisions 


to  leading  order  in  the  parameter  m/>l.  As  will  be  shown  in  this 
paper,  the  results  are  the  same  qualitatively  for  all  impenetrable 
potentials,  but  they  differ  quantitatively  depending  on  whether 
V(x-y)  is  (a)  the  step-function  potential  V'oG :  -x-y  )  with  the 
Heaviside  function  9(x),  or  (b)  the  hard-core  potential  V, wo.  x-y), 
def 1 ned  by 


VHc(x-y)  -rO  for  r>y. 


for  x<y . 


Results  for  any  other  impenetrable  potential  are  reduced  to  those 
in  one  of  these  two  cases.  For  example,  the  soft-core  potential, 
Fsc  (  x-y )  ,  given  by 


Vsc (x-y)  -  e' 


provides  the  same  results  as  VHc(x-y)  .  -  A  penetrable,  localized 
potential  gives  qualitatively  different  results,  but  those  results 
are  the  same  as  those  for  the  delta  potential  where  a  is 
def ined  by 


^  -  j  V (x) dx 


We  will  not  consider  this  case  in  this  paper.) 

In  Sec. II,  we  give  a  variational  expression  for  reflection 
amplitudes  and  relate  it  with  expressions  in  DWBA  and  the  standard 
variational  approach.  We  treat  the  adiabatic  and  anti -adiabatic 


KK 


W5 


regimes  in  Secs. Ill  and  IV,  respectively.  Finally  in  Sec . V ,  we 
discuss  some  possibilities  of  improving  on  the  results  for  the 
cases  of  m  ~  M . 


II.  Variational  Expression  for  Reflection  Amplitudes 


A.  T-matrix 

When  we  divide  the  Hamiltonian  H  into  the  unperturbed  part 
Ho  and  the  interaction  V.  we  can  write  the  Lippman-Schwinger 
equation  for  the  exact  wave  function  j<p,  with  outgoing  wave 
boundary  conditions  as  follows: 


''•>  -  I  gsoi>  +  Go  wv|y',  “  >. 


where  !  *0i>  is  the  normalized  eigen  function  of  Wo  and  Go  is 
the  Green's  function,  defined  by 


1  im  -=—n — : — 

€ -*o*  E-Wo+ie 


Quite  similarly,  the  exact  wave  function  <$}  (_)  I  with  incoming  wave 


boundary  conditions  satisfies  the  following  equation: 


<<po/I  +  <(Pf  ('M  vgq 


As  shown  by  Gross  and  Rungeu,  the  functional  (TD/i  •  introduced  by 


C Tj/,  -  <^o/|vUoi>  +  <^o/[VGo 


-N  A  A  A  A  A  A  A  A  A  A  A  A  />  A\M  VwS'A  A -A\' 


Va’a'a’a' /.'a'a\v\v*a*a' •.*a'a\%\%’a\‘-,.v 


SWa 


+  <£/  H  |  VGo  MV|<p0>> 


-  <Qf  w|VGo  WV-VGo  U)VGo  WVU,  w>, 


is  reduced  to  the  exact  T-matrix  element  <<po/ 1  V  I  ,  when  trial 

functions  |<p,  (^>  and  <y /  |  are  equal  to  the  exact  wave  functions 

|p,(0>  and  |  ,  respectively.  In  addition,  [7]/,  is  stationary 

for  small  variations  of  trial  functions  around  the  exact  wave 
functions.  Therefore  Eq.(9)  gives  a  variational  expression  for 
the  T-matrix. 

When  we  define  |<p,  and  <<pf  ^  |  by 


’>  -  |  <poi>  +  Go  wV|£iW>, 


(10a) 


<<po/l  +  <p/(_)|V,Go  (‘) 


Eq.(9)  can  be  rewritten  in  terms  of  |  <p,  and  <<Pf{~'  I  as 


CHfi  -  <vi  H  |//-£|(p,  (+)>-  <«»o/l//o-£l«>.  <*)>- 


Although  trial  functions  |  <p,  w>  and  <$f  |  are  not  required  to 
satisfy  the  standard  scattering  boundary  conditions,  the  incoming 
(outgoing)  part  of  Ip,  (0>  I  )  should  be  equal  to  !  ;o,'- 

(vO/l).  This  is  the  only  requirement  for  trial  functions  in 

Eq .  ( 1  1  )  . 


ww 
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B.  DWBA 


If  we  know  exact  solutions  of  Ho+V  \  where  V\  is  some 
interaction  potential  appropriate  to  the  problem,  we  can  rewrite 
Eq .  ( 1  1  )  with  the  use  of  the  exact  solutions  of  H  o-V'i,  <yi/  _  ;|  ,  as 
f  ol lows : 


<«>!/  (_)  I  v,  |  <fiQt>  +  <<P({-  \H-E\ot  (*'>  -  <<fi  1/  (*)>. 


where  <v i/  ()|  satisfies 


<<Po/ 1  +  <v  i/  (~  I  V i<7o(* 


The  first  term  in  Eq.(12)  is  the  exact  T-matrix  for  the 
Hamiltonian  Ho+V i . 

When  we  take  \<p,  and  <<p/  as  |  i ,  fo>  and  <<pi/  'I 

respectively,  Eq.(12)  gives  the  DWBA  result: 


C T]/i  -  <<p i/  I  V|  |  <poi>  +  <<P\f  (_)  I  V-Vi  | 


In  this  sense,  the  expression  (12)  (or  equivalently  Eq.vll'  is 
an  extension  of  DWBA. 


C .  Uni  tar  i  ty 

The  variational  expression  for  the  S-matrix  is  given  by  the 
following  matrix  relation: 


C Si  -  1-2-iCT], 


! 15) 


In  general,  however,  this  matrix  [S]  does  not  satisfy  the  unitary 
condition: 

SS*  -  S'S  -  1  .  ('16 ; 

One  way  to  fulfill  Eq.(16)  is  to  formulate  the  problem  in  terms 
of  the  K-matrix13.  but  in  that  case,  we  have  to  use  standing  waves 
i fi  p)  instead  of  q  *'  and  <p(_)  .  One  of  disadvantages  to  use  standing 
waves  is  that  we  have  a  difficulty  in  providing  a  physically 
suitable  form  for  <p  pl  .  This  seems  to  be  very  inconvenient  for 
our  purpose,  because  our  basic  strategy  is  not  to  perform  any 
variational  procedures,  but  to  use  Eq .(11)  to  evaluate  reflection 
amplitudes  with  trial  functions  which  become  exact  in  some  extreme 
limit  of  parameters  involved  in  the  Hamiltonian.  Another 
practica1  difficulty  in  the  k-matrix  formalism  is  that  there 
appears  divergent  integrals  in  the  formula  for  an  impenetrable 
potential . 

Another  way  to  satisfy  Eq.(16)  is  to  ’normalize'  [S]  by  the 
introduction  of  (S] ,  defined  by 

[S]  -  (CSIS*])-,/2CS]  (17a) 

-  C S'j  CS'3CS)r,/:.  176  , 

For  any  matrix  [S] ,  Eqs.(17a  and  17b  give  the  same  fSl .  In  both 


limits  of  m/A/— <=  and  0,  C-S]  coincides  with  C-S]  and  there  is  no 
problem  about  the  unitarity  to  leading  order  in  the  parameter  m  M 


D.  Reflection  amplitudes 

Let  us  apply  Eq .  11)  to  the  problem  of  the  Hamiltonian 

by  Eq .  ( 1  '  )  .  The  normalied  free  state  |po<>  for  the  incident 
is  given  by 

<x  ,  y  |  <p0i>  -  :0)  (x  ,  y  )  , 


with 

0,  ‘0)(x,y)  -  Vh/2  ~k,  exp  ( -i  fc.x)  u,  (  y )  , 
where  L,  is  related  to  the  total  energy  E  through 

E  -  —  +  i  +  - 
C  2m  2’ 

and  Ui(y)  is  the  eigen  function  of  the  harmonic  oscillator, 
by 

u,(y)  -  («/F2‘i  !  )'I/2W,  (y)  exp  (-^y2)  . 


l  ven 

wave 

1  8 ; 


(  19) 


,  20 ) 

l  ven 

21  : 


Here  W,( x)  is  an  Hermite  polynomial.  Similarly,  the  free  state 


With  these  definitions  for  free  states,  the  reflection  amplitude 
C /?]/.  is  given  by 

C /?]/.  -  CSD/, 

-  -2~i  J'dxJdyz/  x.i j  (‘)  (x.y ) 

-2-iJdxJdyz/  •'"  ( x  .  y  1  ( Ho-E)  Vi  ‘(x.y  ('23; 

When  trial  functions  p,  ^(x.y;  =  ■xx.i/l;,  '*J>  and 

so/  (~^*(x.y)  =  <cpf  1 1 ,  y>  satisfy  the  standard  scattering  boundary 
condi t ions 

<p,  ^(x.y)  —  '  ,0)ix.y'  +  X]  #n>0"  c0)(x,y;:  for  x— ®, 

n 

0  for  x-*-®,  (24a; 

and 


<p/  c'‘!*  (x.y)  —  f«J/  (0)  +  ]T]  R/n®n  {0){x,y)  for  x  —  ®, 

i 

0  fori  —  -®  ( 24b 

Eq.(23)  can  be  reduced  to 

CRj/i  -  Rfi-2zi  Jdx  J  dyjf  ( x  .  y ')  ( H-E>  vb  -x.y  .  t 25 ) 

This  is  just  an  expression  in  the  standard  variational 
appr oach .  10  1 1  Therefore  an  extension  of  the  standard  variational 
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i'1-in  ^h.  VittVAVi1 


method  is  made  in  Eq .  (23)  in  which  the  second  integral  replaces 
the  requirement  of  the  conditions  (24a)  and  (24b)  for  trial 
functions. 


III.  Approach  from  Adiabatic  Regime 


In  this  section,  we  evaluate  [/?)/,  in  Eq.(23)  with  trial 
functions  which  become  exact  in  the  limit  of  p>>  1  in  the 
Hamiltonian  (1').  When  p  becomes  very  large,  we  can  neglect  the 
kinetic  energy  term  of  incidon  in  zeroth-order  approximation. 
Then  the  equation  to  determine  the  motion  of  oscillator  can  be 
written  as 


C-i-  +  iy2  +  V  ( x-y )  ]Vn(y  ;  x)  -  Vn(x)Yn(y,x)  , 

<-  e- 


where  the  boundary  conditions  for  the  normalized  wave  function 
Yn( y;x)  are  given  by 


Y„(y\x)  — •  0  as  y  —  ±®. 


In  Eq.(26),  x  is  an  'external'  parameter  and  the  eigenvalue  V„(x) 
provides  and  'adiabatic  potential'  for  the  incidon.  Namely,  the 
motion  of  incidon  is  determined  by 


C-W-  —i  +V'n  (  X  )  )Xn(  X  , 

c-V-  dx~ 


ex„  ( x : 


with  the  following  boundary  conditions: 


-12- 


1 


.•vVv'  -V.V.'V  '  ■  ‘  v  * 


Xn(T) 


A/M/2xfc„e  'k"z 


+  J n/2-kne'k **'0' 


for  r  —  °°  . 


1  0  for  29 

where  S„  is  the  phase  shift  to  be  determined  by  the  solution  of 
the  Schrodinger  equation  (28).  Clearly.  as  and 
Vn(y ; x)  approach  n+1/2  and  un(y)  respectively.  Thus  kn  is 
determined  by  the  relation 


E  - 


2.U 


n  +  2' 


i  30 ) 


By  multiplying  Xn  by  yn.  we  obtain  a  trial  wave  function  as 

<pn  C‘)  ( i ,  y )  -  (fin  H*(x,y)  -  (Pn  (ad,)  -  Xn(x:Yn(  y;x).  (31) 

which  becomes  exact  in  the  limit  y  —  =>  and  is  expected  to  be  still 
a  good  trial  function  for  y>>!  .  The  wave  funcetion  <p„  'ad,J  has  a 
property  of  only  an  elastic  reflection  with  the  amplitude 
exp  ( io„)  . 

The  inelastic  reflection  amplitudes  can  be  evaluated  with 
the  use  of  Eq.(23)  (or  rather  Eq.(25)  in  this  case  ;  and  the  trial 
function  (31).  The  result  is  given  by 

C Rif,  -  e!5'6/, 

-  ^jjp-  / dx Jdy  ( —  (dy ’ 1  ’  V/  ( y ;  * ) C-V,  ‘  (x)X/(x)  +  .V,  v  ;  V/ '  (  x  .  ) 
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^  ’■jc  ^  jc", 


+  ay>(^-'-  ay.^;£),Yt(x)A7(x)}  , 


where  6/,  is  the  Kronecker's  delta  function.  Equation  .32  shows 
that  "non-adiabatic'  transitions  occur  only  when  dV'n  y ;  x  cx  does 
not  vanish. 

In  Fig.l,  we  have  shown  the  numerical  results  of  transition 
probabilities  j  C^?D/o  I J  and  the  normalized  ones  IC&]/olc  for  the  case 
of  the  hard-core  potential  and  p=5.  Solid  curves  represent  the 
exact  results,  given  by  Secrest  and  Johnson15,  while  the  broken 
and  dotted  ones  show  the  results  of  |CR]/ol~  and  |[R]/ol2 
respectively.  We  have  also  calculated  the  case  of  p=0.5,  but  the 
results  for  both  ICR]/il"  and  ICR]/,  |2  are  not  good  except  near  the 
threshold,  namely 


E  ~  max ( i  +  ^ , f  +  ^}  , 


at  which  the  motion  of  the  incidon  is  always  very  slow 
irrespective  of  p. 


IV.  Approach  from  Anti-Adiabatic  Regime 


When  p  is  very  small,  the  kinetic  energy  term  in  the 
Hamiltonian  (1')  becomes  very  important  and  the  free  motion  of 
oscillator  can  be  neglected  in  zeroth-order  approximation. 
However,  we  cannot  neglect  the  potential  term  V  (x-y  ,  because  this 
is  always  larger  than  the  kinetic  energy  for  an  impenetrable 


potential.  Thus  in  zeroth-order  approximation,  we  should  consider 
the  following  equation  first: 


C-^j  +  Xn(xiy)  -  An  ( x ;  y )  , 


with  the  boundary  conditions  for  ,Vn  as 


A„(x;y)  —  f  a/m/2 ~kn  e~'knZ  +  ^/Jx7Z~kne 


for  x-*00, 


for  t  —  -<=. 


It  should  be  noted  that  the  phase  shift  Sn  will  depend  on  y,  in 
contrast  with  the  case  of  Eq.(29).  By  considering  the  difference 
between  Xn(x,y)  and  Xn  ( x  ;  0 )  ,  we  can  write  A„  ( x  ;  y )  in  terms  of 
An  -V  )  ■  Xn(x;  0)  as 


An  (  x  ;  y )  -  e'^Xnix-y)  . 


The  motion  of  oscillator  can  be  determined  in  zeroth-order 
approximation  as 


(-  Jj  ^2  +  2y2^n<'y')  ”  (n+2^'n^y'  • 


This  is  nothing  but  the  equation  for  a  free  harmonic  oscillator 
and  the  solution  is  un( y).  Combining  A„(x;y)  with  un>y),  we  obtain 
a  trial  wave  function  as 


(x,  y)  -  '-fin  (_>*(x.y)  -  <pn'ant,)  (x.  y)  ^un^y  )  e'lk'uXn{  x-y  •  ,  (38) 


IS 


:v_  v\nnn/LTrnrJv.^r.  r  jv.-w 


which  becomes  exact  in  the  limit  y  —  0. 

Since  the  trial  wave  function  (38)  does  not  satisfy  the 
standard  scattering  boundary  conditions  (24),  we  have  to  use 
Eq.(23)  to  calculate  .  The  result  is 

C R3/.  -  ! -2xi J  dx  V (x)Xf(x)Xx(x) 

-  2xi^j-f dx(.X' f(x)+ikfXf(x)2  HXtilde’  ,(x)  +  ik,X,  x)]} 
x  fdy  u/(  y)u,(y)  exp  C-i (ki+k/) y) .  (39) 


In  deriving  Eq.(39),  we  have  used  the  relations 


fdy  u„(y ) u«(y )  e~,au  -  */2‘m!/2nn'.  (-ia)n-me-a2/ilm 


fdy  u„(y)  e~iay  -  i  (|  -  -^)  Jdy  un(y)um,y 


where  L,  a)(x)  is  the  associated  Laguerre  polynomial,  defined  by 


L,  (a)(i) 


_  e  x~°  d* 
m!  dx“ 


(e_Ix~a) 


In  the  following,  we  will  evaluate  Eq.(39)  for  several  types  of 
potentials . 


A.  Hard-core  potential 


When  V (i) -Vhc(x)  ,  Xn(x)  is  easily  solved  to  be 


X„(i) 


![i/ 2-kn  (e  lk'z  -e,k'x)6(x). 


Thus  we  obtain 


[£]/.  ”  -2 Jdy  Uf(y)UiCy)  exp 


C-i  ' 


As  p  approaches  zero  with  the  total  energy  E  fixed,  kn  also 
approaches  zero.  Therefore,  Eq.(44)  has  the  following  expansion 
in  terras  of  p: 


C/?]/i  -  -6/,+2i  */k,kf  I  dy  u/(y)u,(y)y. 


This  agrees  with  the  exact  result  obtained  by  other  methods.02122 
The  elastic  term  becomes  dominant,  while  the  inelastic  terms 
involve  only  one  quantum  jumps  to  leading  order  of  u  whose 
amplitudes  are  proportional  to  . 

In  Fig. 2,  we  have  shown  the  results  for  transition 
probabilities  |(R]/ol2  and  the  normalized  ones  |[A]/ol2  for  the  cases 
of  (a)  m=0.125  and  (b)  m=0.5,  together  with  the  exact  results  of 
Secrest  and  Johnson.15  (Solid,  dotted.  and  broken  curves 
correspond,  respectively,  to  the  exact  results.  l(/?j/ol2>  and 
I  CRJfO  1 2 .  )  The  agreement  between  ICW/ol*  and  the  exact  results  is 
excellent  for  all  calculated  energies  E  for  the  case  of  u=0.125. 
Even  for  m= 0.5,  the  qualitative  features  are  reproduced  in  our 
calculation  of  |(R]/ol2. 
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B.  Soft-core  potential 

We  can  solve  ,Y„(x)  with  the  use  of  the  modified  Bessel 


function  Ku'.z)  for  the  soft-core  potential  Vsc(i).  The  result  is 


-V„ ;  x ) 


'p/2izkn  • 


4kn  (2u/cx2)-'kJa 
a  "  f  ( 1  -  i2fcn/a) 


Kx2k„/a(  ■ 


where  r<  x)  is  the  gamma  function.  We  can  perform  the  integrals 
in  Eq.(39)  and  obtain 


CW]/» 


2 kx  kf 
ki+kf 


(1+M)  O 


a  ' 


xf2£r, (fc.>fc,)/alr(Ui(fc,-Hfc/)/a)  I2  [rd-t-iCfc.-fc/Va)  I2 
la2;  r(l-i2fci/a)  r ( 1  -i  2k//a 'i 

xfdy  uf(y)ui(y)e~,{k''k/)\  (47) 

As  p  approaches  zero,  Eq.(47)  has  the  same  asymptotic 
behavior  as  given  in  Eq.(45).  Namely,  there  is  no  difference 
between  the  hard-  and  soft-core  potentials  in  the  limit  p  — 0. 
The  same  asymptotic  behavior  is  also  expected  for  any  potential 
V7 ( x )  which  becomes  infinite  at  some  point  xo  end  goes  to  zero 
sufficiently  quickly  for  x— <= .  (The  type  of  potentials  considered 
here  is  illustrated  in  Fig. 3.)  This  can  be  seen  rather  easily  by 
evaluating  the  leading  correction  term  in  the  difference  of 
reflection  amplitude,  A(R)/,  ,  between  V(x)  and  Vhc(x-xo)  ■  We  have 


i'9jk  <••*** tAJmXfiti 


, r»4i  < 4 #. » J’ «|| ' 


AC ft]/.  -  -2 rAJ  dxCVCi+xoJ-VHcCi)] 


xc .V/(x)X,(x)  +  ( 1  +p )  e~,fc|IX1  ( x )  £d  t  {  dXJ(t  + 1  k/Xf  ( t )  I  e‘k  ‘ 


+  (\+n)e-,k>IXf(x)J~dt{^±l  +  ifc.A’.CM  }€,m] 


'  J dy  u/  ( y )  u.  ( y )  e*’ 


-i (fc,^k/)y 


where  Xn(t)  is  given  in  Eq.(43).  As  p  goes  to  zero.  Eq .  48;  has 
the  following  form: 


ACf?3/i  -  -4p  <  j  dx  CV(x+xo)  -  Vhc (  x )  ]x 


:  J dy  u/(  y)  u,  ( y )  e~‘  . 


Thus  A[ft]/,  is  at  most  of  the  order  of  p  and  is  negligible. 
Physically,  such  a  general  behavior  stems  from  the  fact  that  as  p 
approaches  zero,  the  de  Broglie  wavelength  of  incidon  becomes  so 
long  that  the  incidon  cannot  see  details  of  the  potential. 

When  we  treat  the  case  of  finite  p  and  small  a.  Eq .  47  will 
give  much  different  results  from  those  of  Eq.(44).  In  order  to 
check  this  situation,  we  have  plotted  |[ft]iol“  in  Fig. 4  for  the 
cases  (a)  a=4  and  (b)  a=0.7  and  compared  with  the  exact  resutls 
of  Secrest  and  Johnson15  and  those  in  DWBA .  vTo  calculate  the 
results  in  DWBA,  we  have  used  Eq . ( 1 4 ■  with  V i ~  -V  x  and 


<P1„  w  (i,y)  -  <p  in  H*  (x  ,  y )  (x  ,  y )  -  u„(y)X„(x),  (50; 

where  An(x)  is  given  in  Eq.(46).)  Solid,  broken,  and  dotted  curves 
correspond,  respectively,  to  the  exact  results,  |(R)ioh  calculated 
with  the  use  of  Eqs.(17)  and  (47),  and  the  results  in  DWBA .  The 
mass  ratio  p  is  taken  to  be  0.5.  For  large  a,  Eq.(47)  gives 
essentially  the  same  results  as  Eq.(44)  and  provides  a  rather  good 
description,  while  DWBA  gives  quite  poor  results.  The  opposite 
is  true  for  small  a.  Physically,  the  heavy  oscillator  will  not 
change  its  motion  appreciably  by  an  impulsive  collision  with  a 
light  particle.  This  is  the  reason  for  the  success  of  Eq.(44) 
for  the  hard-core  potential.  However,  as  the  'interaction  time' 


becomes  longer 

than 

a-'. 

which  can  occur 

for 

the 

very 

sof  t 

potential,  we 

have 

to 

consider 

effects 

on 

the 

motion 

of 

osci 1 lator .  In 

the 

trial 

f  unction 

(38) .  those 

effects  are 

not 

taken  into  account. 

C.  Step-function  potential 

The  step-function  potential  V'o0(-xj-y)  gives  a  little 
different  behavior  for  (R)/,  in  the  limit  p  —  0.  The  wave  function 
Xn(x)  is  given  by 

■Vn(x;  -  a/m? 2~kn  (0(x)  Rne ) 

+  0(-x)  eM}  .  \ 5 1  ) 


-20- 


'y  ^  v v.1*;.' *x*jpij.'v±y  <jv iv 'w-^,-.-.v.v\^vv^. v,v ».- v -<.- 


’'-'WWV’.V 


where 


Rn 


l  kjcHCn 
X  lCn  Kn 


and 


<:  52 


Kn  35  *y2n  (  Vq-E+ti+  1/2)  .  (53; 

As  q  approaches  zero,  [/?]/,  has  the  following  form: 

C R]/l  =  Rid/,  -  ( 1  -/?,)(  1  -rR/) Jdy  u/(y)u,(y)y.  (54) 

Qualitatively.  Eq.(54)  indicates  the  same  results  as  Eq.(45).  In 
addition,  when  Vo  becomes  infinite.  Eq.(54)  is  reduced  toEq.(45). 

The  same  asymptotic  behavior  of  Eq.(54)  is  expected  for  a 
potential  having  the  form  illustrated  in  Fig.5  in  which  the  point 
io  is  determined  by  the  following  relation: 

J  dx  CV(i)  -Vo0(-i+io)I  -  0.  (55' 

The  proof  of  the  above  statement  goes  in  a  similar  way  as  in 
Sec.  IV.  B.  The  difference  A[R]/i  is  given  by  Eq.(48)  in  which  Vhc 
is  replacead  by  VqQ(-x)  and  ,V„  is  defined  in  Eq.(51).  Thus,  as  p 
approaches  zero,  the  leading  correction  terms  for  A[R)/,  is  given 
as  f ol lows  : 

4CR)''  -vm  />  - 


VI  V  K*.  KT 


v0e  ( -i )  ]x 


-  kf(  \-Rf)  (1+R.)  -  fc.ri-R,  l-R/) 


x  /  dy  u/(y  ;u,(y)e' 


■  i (k, ~kf ) y 


In  deriving  Eq.(56),  Eq.(55)  is  used.  As  in  the  previous 
subsection,  ACR]/i  is  at  most  of  the  order  of  p  and  is  negligible 
compared  with  the  terms  in  Eq.(54). 


V.  Discussion 

We  have  applied  a  variational  expression  (Eq.ill))  to 
collisions  of  a  particle  with  a  harmonic  oscillator  in  one 
dimension.  We  have  used  trial  functions  p n  !odi:(i'.yi  Eq.  31;) 
and  '^nfantl)  (i .  y )  (Eq.(38))  in  the  expression  to  obtain  exact  results 
in  the  limits  p->®  and  p  — 0,  where  p  is  the  mass  ratio  m/M.  In 
particular,  we  have  derived  exact  expressions  for  inelastic  as 
well  as  elastic  collisions  to  leading  order  in  the  parameter  p 
for  any  impenetrable  interaction  potential  in  the  limit  p—0. 

For  finite  p,  we  have  obtained  a  rather  accurate  formula 
(Eq.(44))  for  the  reflection  amplitudes,  when  p  is  smaller  than 
unity  and  the  interaction  is  a  hard-core  type  potential.  However, 
for  a  soft-core  potential  and  the  case  of  p~l.  Eq .  i  44  ;  gives 
poorer  results  than  DWBA .  We  can  improve  on  our  results  by 
introducing  some  variational  parameters  in  our  trial  functions 
and  determining  them  by  a  variational  procedure.  One  possible 


way  is  to  give  the  trial  function  in  the  following  way: 


,<  x  .y)  =  l<Pn  (od°(x,y)  +  ;1-A)pn  (ant,)(x.y). 


where  A  is  a  variational  parameter.  This  is  the  combination  of 
the  adiabatic  and  ant i -ad  1  aba 1 1 c  approaches.  Another  way  is  to 
give 


.  i  .  y  =  u„(y )  e',Ak,''J.\nti-Ay  . 


where  \?  x  =  Xn(x'1y=0)  is  defined  in  Eqs.(34:  and  (35^ 


When 


A  =  1  ,  <p,,(x  .  y  ;  is  reduced  to 


(cnt t ,  ■ 


(i,y)  while  for  A=0,  ^n(x,y)  is 


nothing  but  <p„  {DVBA)  (x ,  y)  defined  in  Eq.(50).  Thus  this  is  the 
combination  of  the  ant  1 -adi aba 1 1 c  approach  with  DWBA  . 

In  the  near  future,  we  will  extend  our  method  to  treat 
collisions  in  three  dimension.  in  particular,  gas-surfacde 
scattering  problems. 
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Figure  Captions 


Fig.  1.  Transition  probabilities  for  a  harmonic  oscillator 

interacting  with  a  particle  as  a  function  of  total  energy 
E.  The  hard-core  potential  is  assumed.  The  mass  ratio  p 
is  5.  Solid,  dotted.  and  broken  curves  represent, 
respectively,  the  exact  results  of  Secrest  and  Johnson, 

I  CRD/0 1  “  given  by  Eq.(32'.  and  |[#)/ol2  given  by  Eq.(17)  with 
the  use  of  those  CRj/0- 

Fig.  2.  Transition  probabilities  for  a  harmonic  oscillator 

interacting  with  a  particle  as  a  function  of  total  energy 
E.  The  mass  ratios  p  are  0.125  and  0.5  in  (a;  and  (b), 
respectively.  As  in  Fig.l,  the  hard-core  potential  is 
assumed.  The  exact  results  of  Secrest  and  Johnson  are 
"l,own  by  a  solid  curve,  while  l[R]/ol2  given  by  Eq.(44), 
and  normalized  ones  ICR]/ol'.  given  by  Eq.(17)  with  the  use 

of  those  CRD/o •  are  plotted  by  dotted  and  broken  curves, 

respectively . 

Fig.  3.  A  type  of  potential  which  gives  the  reflection  amplitudes 
in  the  form  of  Eq.(45)  as  p  approaches  zero. 

Fig.  4.  Transition  probabilities  for  the  soft-core  potential  case. 

The  mass  ratio  p  is  taken  to  be  0.5.  Solid,  broken,  and 
dotted  curves  correspond.  respectively,  to  the  exact 
results  of  Secrest  and  Johnson,  ICRJiol2  with  the  use  of 
Eq .  17)  and  (47)  ,  and  those  in  DWBA .  Cases  of  <v-4  and 
a=0.7  are  treated  in  (a  and  b),  respectively. 
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Fig.  5.  A  type  of  potential  which  gives  the  reflection  amplitudes 


in  the  form  of  Eq.(54)  in  the  limit  p— 0. 
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